A mesoscopic description of radiative heat transfer at the nanoscale 
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We present a formulation of the nanoscale radiative heat transfer (RHT) using concepts of meso- 
scopic physics. We introduce the analog of the Sharvin conductance using the quantum of thermal 
conductance. The formalism provides a convenient framework to analyse the physics of RHT at the 
nanoscale. Finally, we propose a RHT experiment in the regime of quantized conductance. 
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It has been discovered in the late sixties that the RHT 
between two metallic parallel plates can be larger than 
predicted using the blackbody radiation form Pn3|. It 
is now known that this anomalous RHT is due to the 
contribution of evanescent waves and becomes signifi- 
cant when the distance separating the interfaces becomes 
smaller than the thermal wavelength Ath = where h 
is Planck's constant, fee is Boltzmann's constant, c is 
the light velocity and T is the temperature. Using the 
framework of fluctuational electrodynamics Polder 
and van Hove (PvH) were able to derive a general form 
of the RHT accounting for the optical properties of the 
media [||. Since this seminal contribution, several re- 
ports have been published in the literature 
quantum-mechanical derivation [12j has confirmed these 
results obtained within the framework of fluctuational 
electrodynamics. While the first papers considered met- 
als, it has been realized that the RHT at the nanoscale 
can be further enhanced for dielectrics due to the con- 
tribution of surface phonon polar itons [H, 14 1. Recent 
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reviews can be found in Rcfs. 

The first attempts to measure a heat flux between 
metallic surfaces at room temperature and micrometric 
distances have proved to be inconclusive [h| lioj . Ex- 
periments in the nanometric regime have clearly demon- 
strated the transfer enhancement 2l|, |22j . Yet the lack 



of good control of the tip geometry did not allow quan- 
titative comparison with theory. More recent experi- 
ments 2^, 2M are performed using silica taking advantage 
of the flux enhancement due to the resonant contribution 
of surface phonon polaritons. A good agreement between 
PvH theory and experiments has been reported 24 1 . 

The purpose of this paper is to establish a link between 
the PvH form of the radiative heat flux and the formal- 
ism of transport in mesoscopic physics. It will help to de- 
velop a more physical understanding of the RHT at the 
nanoscale, which also clarifies how losses and non-local 
effects determine the maximal achievable heat flux (Tol ) . 
Finally, we will show that this reformulation raises the 
prospect of observing quantized conductance for systems 
with sizes on the order of the thermal wavelength Ath- 

We start our discussion with the PvH form of the RHT. 
We consider a vacuum gap with width d separating two 



homogeneous half spaces labeled medium 1 and 2 [see 
Fig. Q] a)]. Then, the heat flux is 0, Eft 



POO J 

Jo 2lT 



)} 



E 



k0 d 2 K (1 



nil 2 



)(1 



n2|2N 



o 4^ |1 
d 2 « 4Im(r^)Im(rf)e- 2Im (T) d 



(1) 



kQ All 2 |l_ r l r 2 e -2Im( T )d|2 



where k = (k x , k y ) and 7 = yfcg — k 2 are the parallel and 
normal wave vector, 0(w, T) = fuo / [cxp(huj / k^T) — 1)] 
is the mean energy of a harmonic oscillator, fco = uj/c, 
rj' 2 arc the usual Fresnel factors for s- or p-polarized 
waves and the sum is carried out over j = s (TE),p (TM) 
and accounts for the two polarizations. In this expres- 
sion, we have clearly separated the propagating (k < fco) 
from the evanescent wave contribution (n > fco). The 
first one can be shown to be described by the usual ra- 
diometric approach 16|. In particular, the term 1 — \r]\ 2 



which is the interface transmittivity is exactly the emis- 
sivity of the interface of medium 1. It is thus tempting 
to extend the definition of the emissivity to evanescent 
waves. This has been proposed in Ref. [13( where Im(rj) 
was introduced as a generalized emissivity. Obviously, 
the generalized emissivity is larger than 1 since it ac- 
counts for an enhanced RHT. A drawback of this ap- 
proach is that it does not provide any physical insight in 
the meaning of an emissivity larger than 1 . 

A different interpretation regarding the evanescent 
contribution only, was given by Pendry Q using 
the single interface local density of states (DOS) 
Im(rj) exp(— 2Im(7)^)| z= o in the space (k,z,u). Here, 
the term local refers to the fact that we introduce the 
DOS in direct space at a distance z above the interface. 
To summarize, the first approach suggests to consider 
that Im(rj' 2 ) is a generalized emissivity whereas the sec- 
ond approach emphasizes the fact that Im(rj' 2 ) is pro- 
portional to the local DOS. 

The interpretation in terms of local DOS is interesting 
as it indicates that the enhancement of the RHT at the 
nanoscale is due to an enhancement of the DOS. Hence, 
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it would be enlighting to develop a formulation of the 
RHT that highlights the number of modes (NOM) in- 
volved. Such a formulation has been developed in the 
context of charge transport in mesoscopic physics and is 
known as Sharvin conductance when the contact area is 
much larger than the electron wavelengths and Landauer 
formalism when the size is on the order of the wavelength 
so that the NOM involved in the transport becomes dis- 
crete. 



AT <C T. Then, we can write the flux <E> as 
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Figure 1: Sketch of the conductance geometry: a) Two ther- 
mal reservoirs are separated by a vacuum gap with width d, 
b) two electron reservoirs with a voltage difference V are con- 
nected through a nanowire. 



Before, reformulating the heat flux in a Landauer-likc 
way, let us remind the structure of the electrical con- 
ductivity in the mesoscopic regime. When considering 
that the dephasing length and the mean free paths are 
much larger than all the relevant length scales, the elec- 
trons can be described by their wave functions. Hence, 
the conduction through a constriction or nanowire [see 
Fig. 1 b)] can be described by a scattering matrix that 
connects the incoming channels to the outgoing channels, 
ft follows that the charge transport can be described by 
means of a sum over all the modes [251 ]2&i 



I=2 i JdEj2 T n(E)[h(E)-f2(E)} 
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where e is the electron charge, T n (E) is the transmission 
probability of a mode n and fi{E) = I/[cxp[(_E — E-p — 
eV)/k B T] + l] and/ 2 (£) = l/[cxp[{E-E F )/k B T} + l] arc 
the Fermi-Dirac distributions and Ep the Fermi energy. 
For small applied voltages V, we obtain 



I = 
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where 2e 2 jh is the quantum of conductance and /o = /2- 
For temperatures k B T E-p, this result reduces to the 
Landauer formula 



(4) 



To derive such a formulation for the RHT case, we 
consider a situation where the two temperatures are close 
enough so that Ti = T + AT and T 2 = T assuming that 
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Here, we have introduced the notation 
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where Dj = 1 — rjrj cxp(2i7<i) is the Fabry-Perot type 
denominator. For propagating modes, Tj 2 {ui, k, d) is the 
two interfaces transmission factor (TF) and takes val- 
ues in the interval [0, 1]. It follows that the heat flux is 
limited by the black-body result. Hence, the expression 
strongly suggests to consider that Tj 2 (uj,K, : d) is also a 
TF for the evanescent modes. As shown in Ref. T 12 
is indeed smaller than 1. This maximal value is reached 
if I'-JH'-Jl = e 2ld which for example can be satisfied for 
s- and p-polarized modes in the case of frustrated total 
internal reflection 27 1 when considering dielectrics and 
also for coupled surface modes (CSM). 

By comparing the expression in Eq. (J3|) with Eq. ([5]) 
one can see that we have already formulated the heat flux 
equation in a very similar manner. The sum over trans- 
verse modes is given by the integral J d 2 K/47r 2 . In order 
to get a manifestly Landauer-like structure, we need to 
integrate the modes over all energies. We first introduce 
the dimensionless variable u = ftw/k B T. Then, by inter- 
changing the order of integration we cast the heat flux in 
the form 
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where n k B Tj3h is the universal quantum of thermal 
conductance [2^,[2^|. Here, we have introduced the mean 
transmission factor (MTF) 



u = f Q °duf(u)T} 2 (u,K,d) 

JT d «/(«) 



(8) 



where f(u) = u 2 e u /(e u - l) 2 and J °° f(u)du = ir 2 /3. 
This is a new quantity which resembles the transmission 
probability for the electrons in the Landauer formula in 
Eq. (|4|). When dealing with electrons, only the TF at the 
Fermi energy is relevant because of the particular form of 
the Fermi-Dirac distribution. By contrast, when dealing 
with bosons, we need to introduce a TF averaged over all 

12 

energies. Let us stress that T- is always smaller than 
1 . This property follows directly from the same property 
of T 12 and the definition of the MTF. Hence, we have 
now a new interpretation of the physical meaning of the 
term involving Im(rj' 2 ) in Eq. ((T|). Instead of interpret- 
ing Im(rj' 2 ) as the local DOS or generalized emissivity of 



Figure 2: Illustration of the contribution of b) propagating 
modes (k < ko), c) modes propagating in the dielectric but 
evanescent in the gap (ko < k < nko) and d) evanescent 
modes confined in the gap (k > nko). 



a single interface, we consider now the two interface sys- 
tem (the gap) and define a TF averaged over all energies 
for a mode with specified n. In this picture, the enhance- 
ment of the heat flux appears as the consequence of the 
increase of the NOM contributing to the RHT. We pro- 
vide in Fig. [2] a) a schematic representation of the modes 
in the (k x ,k v ) plane, from which it becomes clear that 
the number of transverse modes diverges. For n < ko, 
we have the usual thermal radiation due to propagating 
modes, for ko < k < nko, there is a contribution of modes 
which can be viewed as frustrated total internal reflection 
when dealing with dielectrics. Finally, for k > nko, we 
have the contribution of CSM confined in the gap. It re- 
mains to be studied under which conditions the TF takes 
significant values. 



A detailed study of the TF T{ui, k, d) is reported in [31 
The key results are as follows. For a gap thickness d 
larger than the thermal wavelength Ath, the TF is neg- 
ligible for evanescent waves (ko < k) and oscillates be- 
tween and 1 with frequency for propagating waves. It 
is a Fabry-Perot type behaviour. For a gap thickness 
of 100 nm, the TF tends to 1 for both propagating waves 
and waves with frustrated total reflection (ko < k < nko) 
thereby contributing to a significant enhancement of the 
flux. As the TF tends indeed to one for all frequen- 
cies (for s- and p-polarization) , the contribution of these 
modes is simply given by the DOS « 2 ■ (n 2 / 4n) times 
the thermal quantum of conductance in agreement with 
previously reported results (2?1 [32| ■ 

The value of the TF T(uj, k, d) for CSM has a more 
subtle structure. It is essentially negligible everywhere 
in the (uj,k) plane except along the surface mode disper- 
sion relation where it is close to 1 (see Fig. 1 in |31[). 
Hence, after averaging over frequencies, the MTF drops 
by two orders of magnitude as seen in Fig. [3] We also see 
in Fig. H] a) that this drop is smaller for silica than for 
SiC. This is because in the case of silica, there are more 
different surface modes contributing at different frequen- 
cies. Despite the low value of the MTF for surface modes, 



Figure 3: Plot of T p and exp(-2«;d)for 100 nm choosing 
temperature T = 300 K and two SiC [3(| slabs. The modes 
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with k < nko have T p ~ 1. The surface phonon polaritons 
with k > nko have a relatively small MTF but give the dom- 
inant contribution to the heat flux due to the large number 
of contributing modes. The MTF shows for very large k an 
exponential cutoff tx exp(— 2nd). 



a) 1 




Figure 4: Plot of a) T p (d = 100 nm and T = 300 K) and b) 
resulting heat flux $ for SiC and silica [35l ] normalized to the 
black body value $bb = 459.3 Wm 2 . 



the number of additional surface modes is so large that 
the RHT is dominated by the CSM contribution at small 
distance as seen in Fig. [4] b) . 

It is clear from the above discussion that the enhanced 
flux is due to the contribution of surface modes with a 
large value of k. In 31 | 



we discuss the cutoff value of 
k which limits the contribution of the CSM in detail. 
It is roughly given by 1/d as indicated in the graph of 
Fig. [3l A more precise discussion shows [3ll that the 
losses play a key role in defining the exact position of the 
cutoff. It is given for T(u>, K,d) at the surface resonance 
by k > log[2/Im(e)]/d which also sets a cutoff for the 
MTF. 

We finally turn to the ultimate limit of the RHT. So 
far, it seems that the RHT conductance can diverge as 
1/d 2 , since the DOS oc 1/d 2 . Yet, it is known that there is 
a cutoff value for the spatial wave vectors of the phonons 
given by 7r/a where a is the lattice constant. This is 
equivalent to accounting for non-local properties of the 
material HHEl]. 

We now discuss an immediate consequence of our for- 
mulation. If we consider a system as depicted in Fig. [SJ 
we observe that the two wires have a finite transverse 
cross section L x L. Hence, the electromagnetic modes 
are quantized for L on the order of Ath as shown for sur- 
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face plasmon polaritons in Ref. [36| and one expects to 
observe a radiative conductance of the gap with a dis- 
crete NOM similar to the quantized conductance of an 
aperture (37l . l38j |. An experiment similar to the work 
reported in Refs. [39, 40] should be feasible. 
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Figure 5: Sketch of a setup suitable for measuring the quan- 
tized radiative conductance between two wires with a cross 
section L x L. 

To summarize, we have presented a Landaucr-like re- 
formulation for the nanoscale RHT putting it on the same 
footing as the mesoscopic electron transport. In partic- 
ular, we have introduced a mean transmission coefficient 
which is shown to fullfill all the properties needed to give 
a clear understanding of the heat flux on nanoscale in 
terms of the DOS or NOM. In addition, we have proposed 
a near-field heat transfer experiment for measuring the 
quantized conductance. 
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